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Motivation 


• Challenges of simulation-based design 

- High CFD expertise in mesh generation 

► Long setup time 

► High cost due to repeated flow solves on 
fine meshes or high uncertainty due to 
inappropriate meshes 

• Success of error estimation and mesh 
adaptation in goal-oriented simulations 



Objectives 


Adaptive discretization 
of aerodynamic shape optimization problems 


Accuracy 
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Automation 


• Improve design confidence 


• Reduce level of CFD expertise 
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- Direct control over objective 
function discretization error 


- Eliminate the need to handcraft a mesh 
appropriate for all candidate designs 


- Shorten problem setup time 
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Progress toward improved efficiency 

• Reduce cost by systematically increasing depth of refinement 
as designs improve 

- Progressive optimization strategy 
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Problem Formulation 
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min J (X, Q) 
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subject to 


R (X, Q) = 0 
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• Gradient-based optimization 


dJ 

dX 



• Steady Euler equations 


Spatial Discretization: Jh(X,Qh ), Rh(X, Qh) 

• Second-order finite-volume method 

• Cartesian mesh with embedded boundaries 


/ Complex geometry 
/ Automation 
//7-refinement 
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• Gradient-based optimization 


dJ 

dX 



• Steady Euler equations 


Spatial Discretization: Jh{X,Qh), Rh(X,Qh) 


• Second-order finite-volume method 

• Cartesian mesh with embedded boundaries 


/ Complex geometry 
/ Automation 
//7-refinement 




Discretization Error 


Design Space 
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Error Estimate (fixed X) 

J(Q) 
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Number of Cells 


• Leverage adjoint method 

► Error estimates via the method of adjoint weighted residuals 
- Objective function gradient via the discrete adjoint method 
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Dual Role of Adjoints 


Gradients 


H 


Jh = f(X, Q h) 
e.g. C D + (C L -C 


Error Estimates 



e = | Jh ~ Jh) 





Dual Role of Adjoints 



Gradients 


Jh = f(X, Q h) 


d J 


e.g. Cd + (Cl 
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Error Estimation Details 


JhiQh) ~ 


Jc — J h (Q// ) — ll>H ^(Q h) 

Error Estimate 



C Jc-Jh 



Number of Cells 





VH — [^h ~ V* H 


R-/i(Q h ) 


Refinement Indicator 

N 

V = ^2 

1=1 
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Inflow 


Verification: Supersonic Vortex 



Outflow 


Uniform Refinement 





True Error 


Error Estimate 




• No limiter, 0(h 2 ) 

• Effectivity close to 1 
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Inflow 


Verification: Supersonic Vortex 



Uniform Refinement 
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Optimization with Mesh Adaptation 




Modify Geometry 
Mesh 

Flow Solve 
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Evaluate Objective 
Compute Gradient 




Optimize 


• Integration into existing, fixed mesh, 
optimization framework 

- Build sequence of adapted meshes 

- Pass values of objective and gradient 
from finest mesh to optimizer 
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Optimization with Mesh Adaptation 





Modify Geometry 


Initial Mesh 
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Adapt & Solve 


Evaluate Objective 


Compute Gradient 


Optimize 



Integration into existing, fixed mesh, 
optimization framework 

- Build sequence of adapted meshes 

- Pass values of objective and gradient 
from finest mesh to optimizer 


In each design iteration, perform fixed 
(user specified) number of adaptations 

■ Fixed depth strategy 

- Robust and precise control over 
computational resources 
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Optimization with Mesh Adaptation 
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Modify Geometry 


Initial Mesh 
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Adapt & Solve 
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Evaluate Objective 


Compute Gradient 
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Optimize 


• In each design iteration: 

- Start with same initial mesh 

- Adapt until prescribed refinement 
level is attained 

• May be inefficient 



Progressive Optimization 



• Increase mesh refinement in each optimization subproblem 

- Converge a sequence of improving discretizations 


Modify Geometry 


Initial Mesh 


X X* as S 0 


Adapt & Solve 



Evaluate Objective 


Compute Gradient 


Optimize 



15 





Progressive Optimization 


• Increase mesh refinement in each optimization subproblem 

- Converge a sequence of improving discretizations 




Modify Geometry 


X X* as S 0 




Initial Mesh 


Adapt & Solve 



Compute Gradient 



Optimize 





Progressive Optimization 
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• Increase mesh refinement in each optimization subproblem 

- Converge a sequence of improving discretizations 




Modify Geometry 
Initial Mesh 

mmrnmM 

Adapt & Solve 




Evaluate Objective 
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Compute Gradient 


Optimize 
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X X* as £ 


N+2 Cycles 
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Progressive Optimization 


• Increase mesh refinement in each optimization subproblem 

- Converge a sequence of improving discretizations 


• Stopping Criterion 

1 . Gradient or KKT norms, or stall 

2. Specified number of search directions 

3. Diminishing changes in objective 
function 

4. Ratio of design improvement to error: 
refine when T T 

Ji_i - Ji<£ 



Modify Geometry 




Initial Mesh 


Adapt & Solve 




Compute Gradient 


Optimize 



N+2 Cycles 
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Results 


Sonic-Boom Mitigation Inverse Design 

Optimize aircraft shape by prescribing quieter near-field signals 




1 . Pressure-signature analysis 

2. Shape optimization on a fixed mesh 



3. Progressive optimization 
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Pressure Signature of Delta-Wing Body 



Determine pressure signature 3.6 
body-lengths below the model 



Freestream Conditions: 

• Moo = 1 .68 

•Cl =0.15 


h/L = 3.6 
0 = 0 ° 




Mesh and Solution 


Initial Mesh: 12 Adaptations: Isobars 

879 cells 4.5M cells 



Near-field on symmetry plane 


Ap/p 


Pressure Signature 



o—o Initial mesh, 879 cells 
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J (x 0.001) 


Error Convergence 



• Remaining error term is small 
and is 0{h 2 ) 

• Error indicator is 0(h) (due to 
localization) 


Error 


j = 


p 2 

rOC 


J (P-Poo ) 2 dS 


Error bars represent level of 
discretization error 

8 = 2 .J c — Jh\ 



Inverse Design on Fixed Meshes 


Approach: use adaptation to guide construction of a fixed mesh 
for shape optimization runs 4 
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d/dy 


Optimization Targets 





^target) dS 



On-track, 0 = 0° 




Off-track, 0 = 15° 
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d/dy 


Optimization Results 



• 50 design iterations (SNOPT) 

• Ground noise 76.7 PLdB, 9.6 dB reduction in perceived loudness 



On-track, 0 = 0° 




Off-track, 0 = 15° 
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Optimization with Adaptation 

Model Problem Setup 



• Prescribe a target signature from a known shape 

• 1 0 design variables that control body radius 


• Moo = 1.5 and a= 0° 


Initial Shape 
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Optimization with Adaptation 

Model Problem Setup 



• Prescribe a target signature from a known shape 

• 1 0 design variables that control body radius 


• Moo = 1.5 and a= 0° 
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Optimization with Adaptation 


A 



Consider two cases 

1 . Fixed-depth strategy: 7 refinements in each design iteration 

2. Progressive optimization: Increment from 4 to 7 refinements (allow 
designs to advance as far as possible on each level) 



7 Adaptations, ~650k cells 
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Optimization with Adaptation 


Fixed-Depth Strategy Progressive Optimization 

7 Refinements 
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Progressive optimization is about a factor of two faster 
than fixed-depth strategy 


Summary and Outlook 



• Progress toward a gradient- based optimization framework 
with capability to perform adaptive meshing in each design 
iteration 

- Promising approach to enhance accuracy, efficiency and 
automation of simulation-based design 

• Future work 

- Use of error estimates to limit oversolving 

- Transfer of Hessian matrix as the design moves from 
mesh to mesh 

- Dynamic error control and mesh re-use 
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Questions 
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